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We propose a framework to describe and simulate a class of many-body quantum states. We do 
so by considering joint eigenspaces of sets of monomial unitary matrices, called here "M-spaces"; a 
unitary matrix is monomial if precisely one entry per row and column is nonzero. We show that M- 
spaces encompass various important state families, such as all Pauli stabilizer states and codes, the 
AKLT model, Kitaev's (abelian and non-abelian) anyon models, group coset states, W states and 
the locally maximally entanglable states. We furthermore show how basic properties of M-spaces 
can transparently be understood by manipulating their monomial stabilizer groups. In particular we 
derive a unified procedure to construct an eigenbasis of any M-space, yielding an explicit formula for 
each of the eigenstates. We also discuss the computational complexity of M-spaces and show that 
basic problems, such as estimating local expectation values, are NP-hard. Finally we prove that a 
large subclass of M-spaces — containing in particular most of the aforementioned examples — can be 
simulated efficiently classically with a unified method. 



I. INTRODUCTION 

The Pauli stabilizer formalism (PSF) is an important 
tool in quantum information theory. This formalism re- 
gards many-body quantum states, called Pauli stabilizer 
states, that occur as joint eigenstates of sets of com- 
muting Pauli operators. By exploiting the description 
of states in terms of their stabilizers, the PSF provides a 
powerful method to analyze the properties and dynamics 
of stabilizer states in a variety of settings — in fact, the 
PSF is commonly used in virtually all subfields of quan- 
tum information. Important applications include quan- 
tum error-correction [lj, measurement-based computing 
and classical simulations of quantum circuits 3j. In 
addition, the PSF is used in condensed matter physics, 
cf. the study of topological order [4J . 

Notwithstanding its success, a drawback of the PSF is 
that it describes a relatively small class of states. In par- 
ticular, there are only finitely many stabilizer states for 
each given system size. Furthermore these states have 
very particular properties. For example, (modulo some 
trivial cases) they cannot be unique ground states of two- 
body Hamiltonians, every qubit is either maximally en- 
tangled with the rest of the system or completely disen- 
tangled from it, most interesting Pauli stabilizer states 
have zero correlation length etc. [1, H|. In addition, by 
definition the PSF only regards commuting stabilizers op- 
erators. This situation prompts the question of whether 
it is possible to enlarge the class of stabilizer states 
while maintaining a transparent stabilizer-type descrip- 
tion. Such generalizations could lead to new insights in 
many-body quantum states as well as novel applications, 
such as better error-correcting codes, new information- 
theoretic protocols and new quantum states/processes 
that can be simulated efficiently classically. 

A central feature of the PSF is that relevant infor- 
mation about stabilizer states can transparently and effi- 
ciently be extracted by suitably manipulating their stabi- 
lizer groups. Aiming at generalizing the PSF, our goal is 
to identify a mathematical structure which is richer than 



the PSF while maintaining similarly clearcut maps from 
the "stabilizer picture" to the "state picture" . Our route 
towards this end starts with the following observation: 
all Pauli operators are monomial matrices i.e. precisely 
one matrix entry per row and per column is nonzero. The 
basic premise of this work is then to consider arbitrary 
monomial unitary operators (with efficiently computable 
matrix elements) as stabilizer operators, giving rise to a 
general "monomial stabilizer formalism" (MSF). 

Part of the initial motivation for considering the MSF 
stems from the fact that monomial matrices are simple 
operators with favorable mathematical properties. In 
particular, any group generated by a set of monomial 
matrices consists entirely of monomial matrices, making 
such groups rather manageable objects. Most of the mo- 
tivation for studying the MSF, however, is an a posteriori 
one; in fact one of the purposes of the present work is to 
argue that the MSF has several nice features and to make 
a case for the further investigation of this framework as 
an interesting generalization of the PSF. We will do so 
by means of the following two contributions (cf. section 
HIl for a more detailed summary) : 

(a) We show that — perhaps surprisingly — a variety of 
important quantum many-body states are covered by the 
MSF, demonstrating that the latter is significantly richer 
than the PSF. 

(b) We show that basic properties of M-spaces can be 
transparently described by manipulating their monomial 
stabilizer groups; this will lead in particular to efficient 
classical simulations of a subclass of M-states. 

In view of these features, and taking into account the 
wide applicability of the PSF, we believe that the MSF 
provides a promising avenue to describe and simulate in- 
teresting many-body states, potentially leading to new 
applications. Finally we refer to an upcoming work Q 
where the methods developed in the present work arc 
used to arrive at new efficient classical simulations of 
quantum Fourier transforms. 
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II. SUMMARY OF RESULTS 

Next we summarize the main results of this work, out- 
lining in particular the content of contributions (a) and 
(b) mentioned above. 

Generalizing the notion of a Pauli stabilizer code, the 
+1 common eigenspace of a set of monomial unitary op- 
erators will be called an M-space. If an M-space is one- 
dimensional, its (up to a global phase) unique element 
will be called an M-state, generalizing the notion of Pauli 
stabilizer states. 

As for (a) we will demonstrate that the MSF encom- 
passes, in addition to all Pauli stabilizer states and codes, 
the following important state families: the ground level 
of the Affleck-Kennedy-Lieb-Tasaki model @ ; the ground 
levels of Kitaev's quantum double models, which describe 
both Abelian and non- Abelian anyonic systems 0; the 
Laughlin wavefunction at filling fraction v = 1 [9|; the 
family of locally maximally entanglable (LME) states 
[Toj ; coset states of Abelian groups (cf. the Abelian hid- 
den subgroup problem ll l|) : coherent probabilistic com- 
putations [l2|; W states [13I|: Dicke states [14j . 

These examples demonstrate the richness of the MSF. 
They also show that M-states generally do not display 
the aforementioned "special" features of Pauli stabilizer 
states. For example there do exist interesting M-states 
which have non-commuting stabilizer groups and which 
are unique ground states of two-body Hamiltonians. 

As for (b) we will establish basic maps from the stabi- 
lizer picture to the state picture by showing how a des- 
ignated orthonormal basis (called here the orbit basis) of 
any M-space can be constructed when the latter is de- 
scribed in terms of a set of monomial stabilizers. The 
procedure yields an explicit formula for each basis state, 
formulated entirely in terms of manipulations on the sta- 
bilizer group. This result applies to arbitrary M-spaces 
and thus in particular to all examples given above. In 
other words one obtains a single unified method to treat 
a number of seemingly unrelated state families. It will 
also follow from our analysis that all M-states have a 
common, particularly simple structure viz. the nonzero- 
amplitudes of any M-state are all equal in modulus. 

We subsequently use the orbit basis construction to in- 
vestigate classical simulations. Whereas within the PSF 
many quantities of interest can be computed efficiently 
for all Pauli stabilizer states, the situation is different for 
general M-spaces. We will show that one cannot hope 
for general efficient algorithms for several basic problems 
(such as estimating local density operators), as we will 
prove their NP-hardness. In other words the MSF is too 
rich a framework to allow for generally applicable effi- 
cient simulations. However, it is important that these 
results regard worst-case complexity. In fact, based on 
our characterization of the orbit basis, we will identify 
a relevant subclass of M-states for which efficient classi- 
cal simulations can nonetheless be achieved; this subclass 
contains in particular almost all examples listed in (a). 



Throughout the paper we will illustrate our general 
constructions with examples. It is noteworthy that our 
methods allow to recover, with one unified method, the 
classical simulatability of a variety of state families in- 
cluding the Pauli stabilizer states [3J and the quantum 
double models [l5|,[l6j]. In addition, the MSF allows to 
rederive, in a new and unified way, the standard basis ex- 
pansion of stabilizer states in terms of cosets of Z2-linear 
spaces [l7| as well as the matrix product state basis of 
the ground level of the AKLT model Q . 

III. NOTATIONS AND CONVENTIONS 

All Hilbert spaces considered in this work are finite- 
dimensional. We will often consider unnormalized quan- 
tum states in order not to overload the notation. The 
group Q generated by a set of operators U\,..., U m is 
denoted by Q = (Ui, . . . , U m ). We also remark that some 
proofs will be presented in appendices. 

IV. M-STATES AND M-SPACES 

Let T-L be a Hilbert space with orthonormal basis 

B:={\x) :xel}, (1) 

were I denotes some finite set. We will mostly consider 
% to be a multi-party tensor product of d-dimensional 
local spaces and B will usually be a product basis, but 
our arguments hold for arbitrary spaces and bases. A 
unitary operator is called S-monomial if it can be written 
as a product U = PD where D is a diagonal unitary 
operator in the basis B (thus containing phases on its 
diagonal) and where P is a permutation matrix in this 
basis. Equivalently, the matrix representation of U in the 
basis B contains precisely one nonzero entry per row and 
per column. If U and U' are £>-monomial operators, so 
are the operators W and UU' . Furthermore if U is B- 
monomial and V is S'-monomial then U®V is monomial 
relative to the tensor product basis B®B '. 

In the following it will usually be clear from the con- 
text which basis B is considered. Therefore, the prefix 
"£>-" will mostly be omitted i.e. we will simply refer to 
"monomial" operations. Elements of B will typically be 
denoted by |y), \z). 

An important feature of monomial operations is that 
products U1II2 ■ ■ ■ Ut remain monomial regardless of the 
length of the product, as long as every Ui is monomial 
(relative to the same basis). This implies in particular 
that if a group is generated by a set of monomial opera- 
tors, then all elements in this group are monomial. 

The +1 common eigenspace of a set {U\ . . . , U m } of 
operators is the space of all \tp) satisfying 

Ui\ip) — \tp) for every i = 1, . . . ,m. (2) 

A space M. is called a monomial unitary stabilizer space, 
or M-space in short, if there exists a set {Ui} of mono- 
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mial unitary operators with +1 common eigenspace M. 
The group Q — (U-\ , . . ■ ,U m ) is said to be a monomial 
stabilizer group of M. [2a|- A state is called a mono- 
mial unitary stabilizer state, or simply M-state, if there 
exists a set of unitary monomial operators which have 
this state as their unique +1 common eigenvector, up 
to a global phase. In this paper we will only consider 



finite stabilizer groups Q 27|. Otherwise, in the above 



definitions no restrictions are placed on the operators Ui 
except their unitarity and monomiality. 

In the following we envisage M-spaces M. which are 
defined in terms of a set of equations @. A basic prob- 
lem will then be to understand how features of M. can 
be traced back to features of Q. More stringently, we will 
study the computational hardness of determining certain 
quantities (such as expectation values of local observ- 
ables) by suitably manipulating the generators Ui. To 
make meaningful statements about computational effi- 
ciency issues, it needs to be made clear which classical 
descriptions of the Ui are considered to be available (in 
particular one should restrict to M-spaces which have ef- 
ficient descriptions) , and how computational cost is mea- 
sured. This discussion is postponed to section IVIIII In 
sections IVllVIII we will not yet worry about efficiency is- 
sues and our treatment will hold for M-spaces in general. 



V. EXAMPLES 

The paradigm of M-states and -spaces encompasses a 
number of important quantum many-body states: 

• Pauli stabilizer states @ and codes Impor- 
tant examples include the cluster states (l8j . the 
Greenberger-Horne-Zeilinger states [l9j and the 
toric code 

• Generalized Pauli stabilizer states and codes for d- 
level systems [2(| ■ 

• The Affleck-Kennedy-Lieb-Tasaki model Q. 

• Kitaev's quantum double models. These are gener- 
alizations of the toric code which describe systems 
of abelian and non-abelian anyons 

• Laughlin's wavefunction at filling fraction v = 1 

• Locally maximally entanglable (LME) states ficj . 
This class was recently introduced in studies of mul- 
tipartite entanglement. 

• Coherent probabilistic computations [12]. These 
states represent the natural embedding of proba- 
bilistic classical computation (i.e. BPP) into quan- 
tum computation (BQP). 

• Coset states of finite Abelian groups. These states 
occur in important quantum algorithms viz. the 
Abelian hidden subgroup problem (see e.g. 



W-states [l3| and more generally Dicke states [lij ■ 



Remark that the above examples occur in different ar- 
eas of application, ranging from topologically ordered 
systems (cf. quantum double models) to quantum al- 
gorithms (cf. coset states) to multipartite entanglement 
studies (cf. W states and LME states). This richness of 
M-spaces motivates their study as a general framework, 
as initiated in the present paper. We briefly discuss the 
Pauli stabilizer states/codes, the AKLT model and the 
LME states here. See appendix lAl for a discussion of the 
other examples. 



A. Pauli stabilizer states and codes 

An n-qubit Pauli operator has the form a = 701 <E) ■ ■ ■ <E) 
a n where each o~k is either the single-qubit identity oper- 
ator or one of the Pauli matrices and where 7 = ±1, ±i. 
A linear subspace of a n-qubit system is a Pauli stabilizer 
code if there exists a set of commuting Pauli operators 
that has this space as its +1 common eigenspace. If a 
Pauli stabilizer code is one-dimensional then its unique 
(up to a global phase) element is called a Pauli stabilizer 
state. Remark that the Pauli matrices and the identity 
are obviously monomial (relative to the standard basis). 
Since tensor products of monomial matrices are again 
monomial, we find that every Pauli operator is unitary 
and monomial relative to the computational basis. Thus 
every Pauli stabilizer code/state is an M-space/state. 

Remark that general monomial unitary stabilizer 
groups are significantly richer than the Pauli stabilizer 
groups. For example, Pauli operators are product oper- 
ators, whereas elements of general monomial stabilizer 
groups need not be. Second, Pauli stabilizer groups are 
Abelian, whereas general monomial stabilizer groups can 
be non- Abelian (cf. the AKLT model in section IV B[) . 
Finally, the number of Pauli stabilizer groups (and thus 
also the number of Pauli stabilizer states/codes) is finite 
for a given Hilbert space dimension. Monomial unitary 
stabilizer groups on the other hand form a continuous 
family. 

Owing to the richer structure of the allowed stabilizer 
groups, M-states/spaces may exhibit features that can 
never be present in Pauli stabilizer states/codes. For 
example, for every Pauli stabilizer state the bipartite en- 
tanglement between any qubit and the rest of the system 
is either maximal or zero [6j. Also the localizable en- 
tanglement [2lj between any two qubits is either zero or 
maximal 0. It is however easy to give examples of M- 
states where these entanglement measures vary continu- 
ously. It is also known that, except for some uninteresting 
cases, no Pauli stabilizer state can be the unique ground 
state of a two-body Hamiltonian [j| . The example of the 
AKLT model (cf . section IV Bl) shows however that there 
do exist interesting M-states with this property. 
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B. AKLT model 

We show that the ground level of the one-dimensional 
spin-1 AKLT model is an M-space; to our knowledge this 
is a new way of describing the AKLT model. 

Consider two spin-1 particles with local basis 
{|0), |1), |2)}. Let 7r be the projector onto the subspace 
spanned by 



|^> = |01) - |10) 
\rh) = |02) - |20) 



|^) = |12)-|21> 
|V> 4 ) = |00) + |11) 



122) 



(3) 



Now consider a ID chain of spin-1 particles labelled from 
1 to n where one may consider either periodic or open 
boundary conditions. Let -ff^i+i act as I — 7r on spins i 
and i + 1 and remark that -ffj^+i > 0. Let Ai denote the 
zero energy ground state subspace of the AKLT Hamilto- 
nian iJ a kit = HI- Since H^i+x > 0, the space Ai 

coincides with the the common eigenspace of the opera- 
tors Hij+i. In the case of open boundary conditions, Ai 
is four-fold degenerate; in the case of periodic boundary 
conditions, Ai is one-dimensional i.e. the AKLT ground 
state is unique Q. 

Let U be the following monomial unitary operator: 



U 



|01)o-|10) |12) -|21) 

|02) O -|20) 1 00) |11) |22) 1 00). 



(4) 



It is straightforward to show that the +1 eigenspaces of 
U and 7r coincide. Letting L^.i+i act as U on spins i and 
i + 1 it follows that Ai is the +1 common eigenspace of 
the Ui j+i and thus an M-space. 



this immediately implies that \ip) is an M-state relative 
to the product basis B := {U\x)}. 



VI. PRELIMINARY CONCEPTS 

Here we introduce some concepts that will be needed 
in the statement and proofs of our main results below. 



A. The projector p 

Consider a finite unitary group Q with +1 common 
eigenspace Ai and let p denote the orthogonal projector 
onto Ai. Then p coincides with the "group averaging 
operator" 
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Furthermore one has 

Up = p = pU for every U £ Q. 



(6) 



(7) 



For completeness, proofs of these properties are given in 
appendix |Bj 



B. The support of a subspace 



C. LME states 

An n-qubit state is LME (locally maximally en- 
tanglable) if there exists a unitary product operator 
U = Ui ® • • • (§5 U n such that 

C/|V)oc^ 7 *|z) (5) 

for some complex phases 7 X , where the \x) denote n-qubit 
computational basis vectors and where the sum is over 
the entire basis. LME states were recently introduced 
in [l(| where it was shown that \tp) is LME if and only 
if this state can be maximally entangled to an ancillary 
n-qubit system using local controlled-unitary operations 
|28| . It was also shown that the family of LME states 
contains all Pauli stabilizer states. In turn, all LME 
states are M-states. This readily follows from the sta- 
bilizer description of LMEs introduced in [10( , which we 
briefly repeat here. Define the diagonal unitary operator 
D := J2lx\x)(x\ and let Xi be be the Pauli X operator 
acting on qubit i. Since |+) := |0) + |1) is the unique +1 
common eigenvector of the n operators Xi , it follows that 
that \tp') := D\+) n = J2lx\ x ) ^ s the unique +1 common 
eigenvector of the operators Ui := DXiDK Remark that 
the Ui are unitary and monomial relative to the compu- 
tational basis so that is an M-state. As \tp) = U\tp') 



Consider a Hilbert space H with orthonormal basis B. 
The S-support, or simply the support, of a state \^p) is the 
set of basis states \x) E B for which (x\ip) is nonzero. The 
support of a linear subspace Ai of H is the union of the 
supports of its elements. The support of a state/subspace 
will be denoted by supp(|^)) and supp(A / J), respectively. 

Letting p denote the orthogonal projector on the space 
Ai, it is easily verified that the following statements are 
equivalent: 

(a) \x) G supp(X); (b) p\x) + 0; (c) (x\p\x) ? 0. 



C. Uniform superpositions 

Consider a subset O C B and a complex phase 7^ for 
every \y) G O. Consider the states 

W = -m S 7aly) ' I°> = -mE (8) 

V 1^1 \ y )eO Vl^l \ y)e0 

Any state as in the l.h.s. of is called a S-uniform su- 
perposition, or simply a uniform superposition. Remark 
that O is precisely the support of The state \0) is 
called the equal superposition over O. 
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D. The permutation group V 

Consider a unitary monomial operator U — PD where 
P is a permutation and D is diagonal; note that this 
decomposition is unique. Henceforth we will denote U := 
P. Let Q = (Ui, . . . , U m ) be a unitary monomial group. 
The set V = {U : U £ G} is called the permutation group 
of Q. Using that UV = UV for every unitary monomial U 
and V, it can easily be shown that V is a group generated 
by the operators Ui [29j]. The group V naturally acts as 
a permutation group on B. Consider the orbit of |x): 

O x = {\y} : 3P G V s.t. P|ar) = \y)}. (9) 

It is well-known that every \x) belongs to precisely one 
orbit (i.e. the orbit O x ). Straightforwardly applying the 
definition of V one finds that \y) £ O x if and only if 
there exists a complex phase £ and U £ Q such that 
U\x) = £\y). This elementary property will often be used 
in the following. 

Example. Throughout this section as well as section 
IVIII we illustrate the various concepts with a simple "run- 
ning example" . Consider the n-qubit diagonal operator 

T:=a-A® n - A := diag(l,a), a:=e^. (10) 

Further, let Si denote the SWAP gate acting on qubits i 
and i + 1 , where i — 1, . . . , n — 1 . The operators T and Si 
are unitary and monomial (relative to the computational 
basis). We consider the group Q w generated by these 
operators. Since T is diagonal, one has T = I. Since Si 
is a permutation matrix, one has Si — Si. Thus V w is 
generated by the operators Si. It is easily verified that 
V w has n + 1 orbits Oq , . . . , O n , where Oi contains all 
computational basis states \x) where a; is a bit string 
with Hamming weight i (that is, precisely i entries are 
equal to 1). o 

E. The phases £ x (y) 

Let M. be an M-space with stabilizer group Q and let 
p be the orthogonal projector on M. Fix \x) £ supp(A4) 
(hence p\x) is nonzero) and \y) € O x arbitrarily. Then 
there exists U G Q and a complex phase £ such that 
U\x) = We consider the set 5 x ^ y of all phases that 
may occur in this way: 

~ x ^ v :={Z:3UEG s.t. U\x) = £\y)}. (11) 

Remarkably, this set is in fact a singleton. To see this, 
consider £ £ S x ^. y and a corresponding U € Q. Us- 
ing that p — pU it follows that p\x) = £p\y). Thus p\y) 
is proportional to p|x) and the proportionality factor is 
precisely given by £. This shows that this phase does not 
depend on U, so that E x ^ y is indeed a singleton. In the 
following we will denote the unique element of this set by 
£ x (y)- We have proved: 



Lemma 1. For every \x) £ supp(Ai) and \y) £ O x one 
has p\x) = £ x (y)p\y). 

We emphasize that the phase £, x (y) is only defined for 
| a;) £ supp(A^) and \y) £ O x . Finally, note that the 
phase £, x (x) is well-defined for every \x) in the support of 
M, since |x) £ O x . Using that I £ Q and I\x) — \x), we 
in fact find that £ x (x) = 1. 

Example. Consider the group Q w as above and let M w 
be its +1 common eigenspace. We show in section IVIII 
that the basis state |ei) := 1 10 ■ - • 0) belongs to the sup- 
port of Ai w . Note that |ei) belongs to the orbit 0%. The 
other vectors in this orbit are |ea ),..., |e„) where de- 
notes an n-bit string with a 1 in the i-th slot and zeroes 
elsewhere. Letting Pi denote the operator which swaps 
qubits 1 and i (which can easily be obtained as a suitable 
product of Sk gates) one has P % \ei) — |e»). It follows that 
£ei(ei) = 1 for every i. o 

VII. THE ORBIT BASIS 

Consider an M-space Ai specified in terms of a stabi- 
lizer group Q . The latter may e.g. be given in terms of a 
set of generators. Our goal is to construct an orthonor- 
mal basis of M. assuming no prior information about this 
space except for the group Q. In this section we prove 
two results which will directly lead to such a construction. 
The first of these theorems characterizes the support of 
Ai, the second theorem characterizes a designated basis 
of Ai called the orbit basis. 



A. Characterizing the support 

We need some notation. For every \x) £ B let Q x be 
the set of all U £ Q which have \x) as an eigenvector. 
This set is a subgroup of Q. We let {U Xy i, . . . , U Xi i} be 
an arbitrary set of generators [3(| of Q x . 

Theorem 1 (Support of M-space). Consider an Ad- 
space M with stabilizer group Q. There exist orbits Oi 
such that supp(M) = 0\U - ■ -UOd- Furthermore consider 
an arbitrary \x) £ B. Then the following statements are 
equivalent: 

(a) O x C supp(M ). 

(b) (x\U\x) £ {0, 1} for every U £ Q. 

(c) U\x) — \x) for every U £ Q x . 

(d) U x ,i\x) = \x) for every i. 

Proof: lemma[T]shows that for every \x) in the support 
of M and \y) £ O x , one has p\y) ^ so that also \y) £ 
supp(A^). This shows that the entire orbit O x must be 
contained in the support. Thus supp(A'J) is a union of 
orbits as claimed. 
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[a=>c] If O x C supp(vW) then \x) belongs to the sup- 
port of M. since \x) € O x . For every \x) in the support 
the phase £ x (x) is well-defined and in fact £ x (x) = 1 (cf. 
section fVI El) . Consider an arbitrary U € Q x i.e. there 
exists a phase £ such that ?7|a;) = £|x). By uniqueness of 
the phase £, x (x) it follows that £ = £a;(x) = 1. 

[c=^b] Consider an arbitrary U € Q. Then C/|a;) = £\y) 
for some complex phase £ and some If y ^ x then 
(af|Z7|ar) is zero. If y — x then U G Q x and thus £ equals 
+1 by assumption (c). Thus (x\U\x) = £ = 1. 

[b=>a] Let p denote the orthogonal projector onto M. 
Using © we have 



(12) 



Owing to (b) every term (x|Z7|x) in the sum is nonneg- 
ative. Moreover there is at least one nonzero term i.e. 
when U is the identity. This shows that (x|/3|x) is nonzero 
so that \x) belongs to the support of M.. But then O x C 
supp(vW) since the support is a union of orbits. 

[c<S=>d] This equivalence is straightforward since the 
U x .i generate Q x . □ 

Example. Consider the M-space Ai w with stabilizer 
group Gw as in section fVTl For every computational basis 
state \x) one has 



Tlx) 



a ■ a'' 



(13) 



where |x| denotes the Hamming weight of x. This shows 
that, for every x with |x| ^ 1, one has T\x) — \\x) for 
some A ^ 1. Invoking theorem [TJ we find that none of 
the orbits Oi with i ^= 1 are contained in the support of 
A4 W . On the other hand, the orbit 0\ is contained in the 
support. To show this, consider the n-qubit W state 



\W) :-- 



1 



lei) 



(14) 



It is straightforward to verify that T\W) = \W) = S t \W) 
for every i, showing that \W) € M w . Since this state 
has support 0%, it follows that 0\ C supp(A4) K1 . In 
conclusion, the support of M. w is identical to 0\. o 



Lemma 2 (Orbit states). Every \ip x ) is a uniform su- 
perposition state given by: 

\^ = T^TT E Z*(V)\V)- (16) 
Proof: using that p is given by (|6|) one has 

lVx> c< J2 u \ x ) = E c v\y^ 

ueg |y)eOx 



for some coefficients c y . This shows that supp (!?/>?)) Q 
O x . Furthermore for every \y) € O x there exists U € Q 
such that U\x) = £, x (y)\y). Using that U'\ip x ) = \ip x ) it 
follows that 



(v\ipx) = £.x(y) ■ (x\i/j x ). 
In combination with (|17l) this implies that 

\i/) x ) = (x\ip x ) &(y)\y)- 

Since \ip x ) is normalized, it follows that 

a 



(x\ip x 



(18) 



(19) 



(20) 



for some complex phase a. Finally, it follows from def- 
inition (fT5j) that (x\%[) x ) = \\p\x)\\ > 0. This shows that 
a=l. □ 

Since \tp x ) is a uniform superposition over an orbit of 
V, we call this state an orbit state. More precisely we 
call \ip x ) the orbit state determined by \x). 

Example. The orbit state of M w determined by |ei) is 
the n-qubit W-state: 



E MiOlv) 




= \W). 



(21) 



B. Constructing the orbit basis 

For every \x) G supp(.A/f) consider its normalized pro- 
jection onto M: 



p\ x ) 

\W)\\ 



(15) 



Henceforth whenever considering a state (fl"5)) we will tac- 
itly assume that |x) belongs to the support since other- 
wise p\x) = 0. Note that by construction U\ip x ) = \ij) x ) 
for every U G Q. Interestingly, the states \ip x ) have the 
following explicit form: 



Next we show that a basis of M. can be constructed 
by selecting a suitable subset of orbit states; this is the 
orbit basis. 

Theorem 2 (Orbit basis). Consider an M-space A4 
with stabilizer group Q and support supp(M) = 0\ U 
• • • U Od- Choose an arbitrary representative \xi) in each 
orbit Oi and the associated orbit state \ipi) :— \ip Xi ) for 
every i from 1 to d. Then the set ^ := • • • , IV'd)} 

is an orthonormal basis of A4 . Furthermore this basis is 
independent (up to global phases) of the choice of orbit 
representatives \xi). We call \& the orbit basis of M.. 
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Proof: since B is a basis of the Hilbert space and since 
\ip x ) is defined as the projection of |x) onto M., the col- 
lection of all orbit states span M. (although generally 
these states are not linearly independent). Consider an 
arbitrary |x) in the support of M.. By construction there 
exists an i between 1 and d such that |x) £ ft. Lemma 
[T] shows that p\x) oc p\xi) so that \ip x ) c< IV^)- Since the 
orbit states span A4, it follows that the states \ipi) span 
M. as well. Since \ipi) has the orbit Oi as its support ow- 
ing to (|16p and since these orbits are mutually disjoint, 
the states \ipi) are mutually orthogonal. This shows that 
^ is an orthonormal basis. Finally, lemma [T] shows that 
the basis $ is independent (up to global phases) of the 
choice of orbit representatives |xj). □ 

Theorems [T] and [2] show that the following procedure 
allows to correctly identify the orbit basis of A4 . 

• Determine all orbits of V and consider a represen- 
tative Xfc in each orbit Ok- 

• For each k, decide whether Ok Q supp(A4) by 
means of the characterization in theorem [TJ 

• The orbit basis ^ is the collection of all orbit states 
\tpxk) f° r which Ok C supp(A^). 

As desired, this procedure can be implemented by means 
of manipulations on the stabilizer group Q . Furthermore 
each orbit state is itself characterized completely in terms 
of properties of Q owing to lemma [5J 

Example. We have seen that supp(A / f UI ) = 0\. The- 
orem [5] thus shows that M. w is one-dimensional, with 
|^ei) as its unique element. Furthermore we have shown 
that \ip e ±) = \W). Thus, the W state is an M-state with 
stabilizer group Q w . o 

Remark that the orbit basis construction applies to 
arbitrary M-spaces and thus in particular to all instances 
given in section |Vl This yields one unified method to 
analyze all these state families; see section HXl for further 
illustrations. 



C. Some corollaries 

Next we discuss some immediate corollaries of theo- 
rems [T] and [2j The first corollary is interesting in that 
bounds the dimension of M. by means of a purely combi- 
natorial quantity viz. the number of orbits. Remark that 
there exist well-developed tools to count /estimate orbits 
of permutation groups, which may thus be imported into 
the study of M-spaces. 

Corollary 1 (Dimension). The number of orbits d con- 
tained in the support of M. (cf. theorem^) coincides with 
the dimension of M.. It follows that this dimension is up- 
per bounded by the total number of orbits of V . 



The above corollary will be used in section IIXI to give 
a simple proof that the AKLT ground level in the open 
boundary conditions case is four-fold degenerate. 

Second, theorem [2] leads to the following characteriza- 
tion of M-states. 

Corollary 2 (M-states). Every M-state \tp) is a uni- 
form superposition. More precisely, if |x) is an arbitrary 
basis state satisfying (x\ip) ^ 0, then l^) oc \ip x ) where 
\ip x ) is given explicitly in lemma\^ 

Remark that the support of any M-state coincides with 
precisely one orbit. This implies in particular that knowl- 
edge of a single |x) satisfying (x| , 0) ^ implies complete 
knowledge of the entire support of which then must 
coincide with O x . Theorem [T] may be used to determine 
a suitable x such that 1-0) oc \ip x )- 

It is interesting to compare corollary [2] with a charac- 
terization of Pauli stabilizer states obtained in [l7[ • Con- 
sider a Pauli stabilizer state on n qubits. Then there 
exists a linear subspace S of Z£ , an x € Z£ and complex 
phases £j, such that 

m = -L E (22) 

v ' ' yes+x 

Thus every Pauli stabilizer states is a uniform superposi- 
tion, the support of which is identified with a coset x + 5. 
Corollary [2] shows that, in fact, all M-states have a ba- 
sis expansion with an analogous structure. Corollary [2] 
can in fact be used to rederive (|22[) in a simple way; see 
section IIXI 

Restricting attention to many-qubit systems there is 
an noteworthy connection to LME states. In the termi- 
nology of the present work, an n-qubit state is LME 
iff there exists a product basis B such that this state is 
a S-uniform superposition where suppdV 1 )) is the entire 
basis B. Corollary [2] shows that n-qubit M-states (when 
monomiality is considered relative to product bases) can 
be regarded as generalizations of LME states where uni- 
form superpositions with arbitrary supports are consid- 
ered. 

Finally, we discuss a subclass of stabilizer groups for 
which the results above can be simplified. A monomial 
unitary group Q is called pure if it has a generating set 
of the form {Pi, . . . ,Pk 7 Ai, . . . , A m } where every Pi is a 
permutation and every Aj is diagonal. Such a generating 
set is also called pure. For every |x) let \O x ) denote 
the equal superposition over the orbit O x (recall section 

rvrci) . 

Corollary 3 (Pure stabilizer groups). Let A4 be an 

M-space with pure stabilizer group Q . Then the orbit basis 
ofM has the form * = {|0 X1 >, . . . , \0 Xd )}. 

Proof: let V be the permutation group of Q. Consider 
a pure generating set {Pi, Aj} of Q. Since Pi — Pi and 
Aj = I, it follows that V is generated by the operators 
Pi. This implies that V C Q. Let |x) belong to the 
support of A4 and consider an arbitrary | y) G O x . Then 
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there exists P £ V such that P\x) = \y). Since P e Q 
and owing to the uniqueness of the coefficient £, x (y) it 
follows that £ x (y) = 1- Invoking lemma [2] it follows that 
\rp x ) = \O x ). □ 

Examples of pure stabilizer groups are numerous: con- 
sider e.g. the toric code states, quantum double models, 
W-states, coherent probabilistic computations and coset 
states of Abelian groups (cf. section [V] and appendix 
|A|) . We will consider the example of quantum doubles in 
more detail in section HXl 



VIII. COMPUTATIONAL COMPLEXITY AND 
CLASSICAL SIMULATIONS 

Thus far, we have studied general mathematical fea- 
tures of M-spaces, not worrying about the computational 
complexity of the concepts involved. Here we discuss 
such issues. For concreteness, we consider many-qubit 
systems where monomiality is defined relative to the 
computational basis {|x)}; generalizations are straight- 
forward. 



A. Classical descriptions of M-spaces 

First it needs to be made clear which classical descrip- 
tions of M-spaces are considered to be available. It is 
natural to consider n-qubit M-spaces M. that meet the 
following requirements: 

(i) M has a stabilizer group Q = (Ui, . . . , U m ) with 
m — poly(n) generators. A classical description of 
each Ui is considered to be given as an input; each 
of these descriptions is assumed to be efficient. 

(ii) Every generator Ui has efficiently computable ma- 
trix elements in the following sense. Suppose that 
Ui acts on the computational basis as 



Ui 



\i(x)\lTi{x)} 



(23) 



for some complex phases Aj(x) and for some per- 
mutation TTi of the set of n-bit strings. Then Ui is 
said to be efficiently computable if, given any n-bit 
string x as an input, the following two conditions 



arc fulfilled 31 



— There exists a poly(rt) time classical algorithm 
to compute 7Tj(x) and 7r~ 1 (x). 

— There exists a poly(n, k) time classical algo- 
rithm to compute the phase (x) up to k bits 
of precision. 

These conditions entail that, given any row of the 
matrix Ui , it is possible to efficiently compute which 
matrix element within that row is nonzero and what 
the value of that matrix element is, and a a similar 
condition for the columns. 



The above conditions are met in many cases of interest. 
In particular we have: 

Except for the LME states, all M-spaces considered in 
section\V\ satisfy (i)-(ii). 

This statement is easily verified and the arguments are 
omitted here. LME states generally do not satisfy (i)-(ii) 
since the phases ^ x in ([5]) may not be efficiently com- 
putable. However, all LME states where the function 
x — > 7 X is classically computable in poly(n, k) time up 
to k bits — which is a large class — do satisfy conditions 
(i)-(ii). 



B. Computational complexity 

Considering n-qubit M-spaccs described as above one 
may investigate the computational hardness of a variety 
of tasks; here an algorithm is considered to be efficient if 
it runs in poly(n) time. Natural problems are: 

(PI) Decide if the +1 common eigenspace A4 of the gen- 
erators Ui is nontrivial. 

(P2) Given an M-state sample classically from the 
distribution {Kyl^)! 2 }- 

(P3) Given an M-state, compute the expectation value 
of a fc-qubit observable for some constant k. 

We will show that one cannot hope for efficient classical 
algorithms for PI-P3 that apply to all M-spaces in gen- 
eral. This is a point where the MSF sharply contrasts 
with the Pauli stabilizer formalism, where many prob- 
lems of interest can be answered efficiently for arbitrary 
Pauli stabilizer states and codes; this holds in particu- 
lar problems PI to P3. Hardness of the above problems 
shows that one should look for relevant subclasses of M- 
spaces for which efficient solutions are possible. 

The intractability of P1-P3 in fact holds even for very 
simple M-spaces viz. those with diagonal, local genera- 
tors Ui. 

Problem 1. The input is a set of n-qubit diagonal uni- 
tary operators {U\, . . . , U m }. Each Ui acts nontrivially 
on at most 3 qubits and all matrix entries of Ui (in the 
computational basis) are either 0, 1 or —1. The problem 
is to decide whether these operators have a +1 common 
eigenvector. 

Problem 2. The input is a set of diagonal unitary 
operators Ui as in Problem 1, with the following addi- 
tional constraint: it is promised that these operators have 
a unique (up to a global phase) +1 common eigenstate 
which is thus an M-state. The problem is to sample 
classically from the distribution {|(x| , i/')| 2 }- 

Problem 3. The input is as in problem 2. The prob- 
lem is to compute (ip\Zi\ip) with accuracy e = l/poly(n), 
where Zi is the Pauli a z operator acting on qubit i. 
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Theorem 3. None of the problems 1-3 can be solved 
classically in polynomial time unless P = NP. 

Theorem [3] will be proved by reductions to (variants 
of) the satisfiability problem, which is NP-complete. See 
appendix [Cj 



C. Classical simulations 

We focus in more detail on the classical simulation 
problems P2 and P3. Even though these tasks are in- 
tractable in their worst case, efficient solutions exist for 
subclasses of M-states. Here we provide sufficient criteria 
for the existence of efficient algorithms. 

Theorem 4. Consider an n-qubit M-state \tp) . The sta- 
bilizer group Q is described in terms of m= poly(n) ef- 
ficiently computable generators {U\, . . . , U m } as above. 
Let A be a k-local observable with k = O(logn) and 
\\A\\ < 1. Suppose that the following tasks have efficient 
classical algorithms: 

(a) Determine any \x) such that \ip) oc \il> x ); 

(b) Generate a uniformly random element in O x . 

(c) Given \y), decide if \y) G O x . 

(d) Given \y) G O x , compute £, x (y). 

Then there exists an efficient classical algorithm to sam- 
ple the distribution II := {|(y|^)| 2 }. Furthermore, then 
there exists a efficient classical algorithm to estimate 
(ijj\A\ip) with accuracy e = \/poly{n) with success proba- 
bility that is exponentially close to 1. 

Proof: since is a uniform superposition with sup- 
port O x for some \x) (recall corollary [5]), II is the uniform 
distribution over this orbit. This shows that efficient clas- 
sical algorithms for (a) and (b) imply an efficient algo- 
rithm to sample II. 

We consider the second claim. Every fc-local observable 
can be written as a linear combination of poly(n) Pauli 
operators where each coefficient in the linear combination 
has modulus not greater 1. Therefore it suffices to prove 
the claim for Pauli operators . For every a G ZJ? consider 



X(a) 
Z{a) 



X ai 



)X a„ 

z a ~. 



(24) 



Then every Pauli operator can be written as r yX(a)Z(b) 
for some a, b and some 7 € {±l,±i}. We refer to ap- 
pendix [Elfor some standard properties of Pauli operators. 
For every |y) € O x define 



t [V) ~ I otherwise. 1 > 



Using that \ip) oc \ip x ) and lemma [2] one finds 



(26) 



By a standard Chernoff bound argument (see e.g. the ap- 
pendix of |23j), the sum in (|2^|) may be estimated with 
l/poly(n) error with exponentially small failure proba- 
bility by generating poly(rt) random elements \yi) G O x 
and by computing the average of F(yi). Furthermore 
owing to assumptions (a)-(d) this procedure can be im- 
plemented in polynomial time. □ 

Remark that (a) can be approached with theorem [1] 
Furthermore, in order to compute £ x (y), it suffices to 
determine any single U € Q such that U\x) oc \y), since 
then £, x (y) is simply given by the matrix element (y\U\x) 
(cf. section tYl Ej) . Finally, we point out the following 
elementary approach to (b): 

Lemma 3. The following procedure generates a random 
\y) G O x . First, generate a random permutation P G V . 
Then compute \y) = P\x) and output y. 

The proof of lemma [3] uses basic group theory argu- 
ments and is given in appendix |DJ Interestingly, there 
exists a well-developed theory of (approximately) gener- 
ating random elements in finite groups and efficient al- 
gorithms are available for a variety of groups (see e.g. 
|2 ij ) . These methods may be thus imported to the study 
of classical simulations of M-states. 

Efficient algorithms for (a)-(d) exist for a variety of 
M-states. In fact it can be shown that: 

Except for the LME states, all M-states considered in 
section\V\ satisfy (a)-(d) in theorem^ 

As before, LME states for which the function x — » 7 X 
can be computed efficiently do satisfy (a)-(d). In sec- 
tion |IX] we work out the examples of stabilizer states and 
quantum double models. 



IX. APPLICATIONS 

Throughout this paper we have illustrated our results 
by means of the simple example of the W states. Here 
we give some more sophisticated examples viz. the Pauli 
stabilizer states, the AKLT model and Kitaev's quantum 
double models. 



A. Pauli stabilizer states 

We show how the MSF can be used to rederive some 
interesting features of Pauli stabilizer states in a new way. 
In particular: 

(i) We rederive the expansion (|2"2"|) first proved in [ItJ ■ 

(ii) We show that the conditions (a)-(d) of theorem 
|4] are fulfilled, thus showing that Pauli stabilizer 
states can be efficiently simulated classically in the 
sense of theorem |4j this recovers (a variant of) the 
result Q. 
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Let |V) be an n-qubit Pauli stabilizer state with Pauli 
stabilizer group Q . Recall that every minimal set of gen- 
erators of Q contains precisely n Pauli operators [22j . 
We arbitrarily fix such generators {<7i, . . . , <7 n }. Since 
Y = iXZ, every generator can be written as 0$ = 
'jiX (s l ) Z (t l ) where ji G {±l,±z}; recall also the no- 
tation (|24l) . Remark that o~i is a product of a permuta- 
tion matrix X(s l ) and a diagonal matrix jiZ(t l ) so that 
<7j = X(s l ). Let S be the Z2-linear space generated by 
the vectors s\ Using (|E1IE3|) one easily finds that 

V = {X(s) : s e S}; O x = {\x + s) : s G S} (27) 

for every i £ ZJ. Invoking corollary [5] these identities 
immediately imply that 1-0) has the form (|22l) for some 

We address (ii). First, we show that an x such that 
\ip) oc \ifj x ) can be computed efficiently. Recall the defini- 
tion of the group Q x given in section [VII Al Let V denote 
the subgroup consisting of all a G Q satisfying a cx Z(b) 
for some b. It follows from (HU [E2]) that Q x = V for 
every x. We now claim: 

Lemma 4. A generating set {D\, ■ ■ ■ ,D{\ of T> can be 
determined in poly-time; moreover Dj = (— l) Uj Z(d J ) for 
some (efficiently computable) Uj G Z2 and d? € Z2 ■ 

This result can be proved using standard Pauli stabilizer 
arguments; for completeness a proof is given in appendix 
lEl Theorem [T] now implies that O x — supp(|V)) if and 
only if Di\x) = \x) for every i. Using (|E2|) this is equiva- 
lent to requiring that x d? = Uj for every j. A solution 
x to this system of equations can be computed efficiently. 

Since we have access to a generating set {s 1 , . . . , s"} of 
S, we can efficiently determine whether y G S + x, given 
y as input. Also a random element in x + S can easily 
be generated efficiently. 

Finally we show that, given any y G x + S, the phase 
£ x (y) can be computed efficiently. First we compute an 
arbitrary a — (<2i, . . . , a n ) G Z?? satisfying ^ a,;s' = x + 
y; this regards solving a system of linear equations and 
can be done efficiently. Using the properties (|EHIE5[) it 
follows that for such a one has 

at 1 ...<" cx X£> S *)Z(J>f) 

(28) 

= X(x + y)Z(j:a i t i ). 

It follows that cr" 1 ...o-^ n \x) = £|y) for some complex 
phase £; by the uniqueness of £ x (y) w e have £ = £,x(y)- 
Given a, x and ?/ it is straightforward to compute £ in 
poly-time. 

B. AKLT model with open BCs 

Consider the AKLT model with open boundary condi- 
tions. Let A^ opcn denote the ground level subspace. We 
will use the MSF to prove the following properties, first 
proved in [8|: 



(i) The ground level is 4-fold degenerate; 

(ii) An orthonormal basis of ground states is given by: 

IV'opcn) = Tr{crer ai . . . cr Q „ } |o x . . . a„). (29) 

Here a G {7, X, Y, Z}, the sum ranges over all a% G 
{0, 1, 2} and <7o := X, a\ := Y and 02 : = Z. In fact 
we will prove that ([2^]) is the orbit basis of A4 open . 

We prove the results for n even. Let g°P cn be the 
group generated by the operators t/^+i. Let P denote 
the permutation operator obtained by replacing all minus 
signs in (j4j by plus signs. Then obviously U%,%+i = P%,i+i 
so that p°P° n is generated by the operators Pj,j+i. If n 
is even then it is straightforward to show that p°P cu has 
4 orbits: 

• Oi contains all basis states with an even number 
of Os, an even number of Is and an even number of 
2s. A representative element is \aj) = |0 . . . 000). 

• Ox contains all basis states with an even number 
of 0s, an odd number of Is and an odd number of 
2s. A representative element is \ax) = |0 . . . 012). 

• Oy contains all basis states with an even number 
of Is, an odd number of 0s and an odd number of 
2s. A representative element is \ay) = |1 • • • 102). 

• Oz contains all basis states with an even number 
of 2s, an odd number of 0s and an odd number of 
Is. A representative element is \az) — |2 . . . 201). 

The Pauli matrices satisfy the commutation relations 

OWi = — (TiGo ClCa = — O"201 /oq\ 

OW2 = -0200 0000 = 0101 = 0202(= I)- 

Using these relations it readily follows that each IV'open) 
is a +1 eigenstate of all operators f/i,i+i, so that these 4 
states belong to A^ open . Using the commutation relations 
([50]) and the fact that X, Y and Z are traceless operators, 
one also finds that IV'open) nas support O a . This shows 
that these four states are orthogonal. Therefore the di- 
mension of A'! ' 3611 is at least 4 and each orbit 0„ belongs 
to the support of A4 opcn . Owing to corollary Q] the di- 
mension is at most the total number of orbits of "P°P en ; 
which equals 4. This shows that .A/f opon is 4-dimensional 
with basis (j2"9")l . 

Finally we show that IV'open) i s the orbit state deter- 
mined by \a a ). Since all four orbits O a belong to the 
support of .M opon , the orbit basis consists of the four 
corresponding orbit states. As IV'open) belongs to J\4 open , 
this state must be a linear combination of the states in 
the orbit basis. Since |Vo P en) nas support O a as argued 
above, the latter is only possible if IV'open) * s the orbit 
state determined by \a a ). 

By using the commutation relations (|30[) and by ap- 
plying the definition of the phases £ x {y) it can also be 
shown directly that |Vopen) is the orbit state determined 
by \a a ). This argument is omitted here. 
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C. AKLT model with periodic BCs 

Consider the AKLT model with periodic boundary 
conditions and n even. Let A4 pcl denote the ground 
level subspace. We rederive the following properties, first 
proved in 0: 

(i) The ground level is non-degenerate; 

(ii) The unique ground state is (using an analogous no- 
tation as in (|29|l ) 

IVpcr) oc ^Tr{cr ai . ,.o- 0n }|ai • ..On). (31) 

Using the commutation relations (|30l) it is straightfor- 
ward to show that \ip pe i) belongs to M pcr . As above, 
•pper j g g enera ted by the permutation matrices Pi.i+i 
(now imposing periodic boundary conditions). Further- 
more one may verify that V pcr has the same four orbits 
as in the open boundary conditions case. Note also that 
|V>per) has support Oi so that 0/ C supp(A4 per ). We 
show that none of the other orbits are contained in the 
support. To do so define the operators 

A = U 

B = U n - 2 ,n-iA. (32) 

One can then verify that B\a a ) = — \a a ) for every 
a = X,Y, Z. Since B g G pm , owing to theorem Q] this im- 
plies that none of the orbits Ox,Oy or Oz are contained 
in the support of Ai peT . It follows that the support coin- 
cides with the single orbit Oj. Invoking corollary [1] then 
implies that Ai pcr is one-dimensional. 

D. Quantum double models 

Consider a quantum double model defined on a sphere 
(cf. appendix IA2[) . Such a system has a unique ground 
state IV'qd) \Mi which is thus an M-state. Letting Q be 
the stabilizer group of this state, here we show: 

• Q meets requirements (a)-(d) of theorem[U showing 
that |^qd) can be simulated classically in the sense 
of theorem [4] 

The reference work [|| can also be used to show that 
standard basis measurements can be simulated classically 
(although it is not explicitly discussed there). Efficient 
simulations of local expectation values were previously 
achieved using tensor network methods [IH fl6| . 

First we describe the stabilizer group. Recall the defi- 
nition of the operators U p and V v (k) defined in appendix 
IA 21 Since U p is diagonal and V v (k) is a permutation ma- 
trix, the group Q generated by these operators is pure. 
Note that the U p mutually commute since these are di- 
agonal operators. Furthermore it is easily verified that 
[V v (k),U p ] = and [V v (k),V w (l)] = for all vertices 
v w, for every plaquette p and for every k 7 l G G. 



It follows that a general element U G Q has the form 
U = PD where 

P = T\V v (k v ) where k v G G 

V 

D = YI U p" where x p G {0, 1}. (33) 
p 

Since P is a permutation matrix and D is diagonal, one 
has U = P. 

Next we determine the support of IV'qd}- A standard 
basis state has the form |g) = \g e ) where the product 
is over all edges e and where g e G G. Let S be the set 
of all |g) satisfying U p \g) = |g) for all plaquettes p. We 
claim that supp(lV'qd)) = S. To see this, first note that 
theorem [T] shows that supp(|-0 q d)) ^ S. To prove the 
reverse inclusion, consider an arbitrary U = PD G Q as 
in (J53J). If |g) € «S then D|g) = |g). Consequently, 

<g|tf |g> - (S\PD\ S ) = (g|P|g) G {0, 1}, (34) 

where the last inclusion holds since P is a permutation 
matrix. Invoking theorem [T] proves the claim. 

Recall that |"0qd) is an M-state and that Q is pure. 
Corollary [3] implies that l^qd) is the equal superposition 
over 5, and that 5 must be an orbit of V . Since one 
manifestly has |e) G S (where e has the neutral element 
e in all its entries), it follows that S — O e . Consequently, 
IV'qd) coincides with the orbit state IVO- This shows that 
(a) in theorem 0] is fulfilled. 

To show (b) we use lemma[3J Since an arbitrary P G V 
has the form given in (|33p . a random P can efficiently 
be generated by generating a random k v G G for every 
vertex v. Applying P to |e) yields a random state in O e . 

Condition (c) holds since the support O e = S is defined 
in terms of poly normally many constraints U p \g) = |g), 
each of which is easily verified. Finally, since IV'qd) is an 
equal superposition state, one has £ e (g) = 1 for every 

g e O e . 
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Appendix A: Examples of M-states and -spaces 

1. Pauli stabilizer formalism for qudits 

Generalizations of Pauli X and Z matrices exist for 
d- level systems ("qudits") with arbitrary d as follows: 

X d = £> + l)<ar|; Z d = £ e 2 ™ ' d \x){x\. (Al) 

x=0 x=0 
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Here {|0), . . . \d — 1)} is the standard basis and x + 1 
is computed modulo d. Based on these definitions one 
may consider generalized Pauli operators acting on n qu- 
dits and the associated stabilizer states and codes. It 
is straightforward that generalized Pauli operators are 
also monomial unitary operators relative to the standard 
basis, so that stabilizer states and codes for qudits are 
M-states and-spaces, respectively. 

2. Quantum double models 

Consider a finite group G and a Hilbert space "Hioc 
with basis {\g) : g £ G} where basis vectors are labeled 
by group elements. Furthermore consider a 2D square 
lattice A arranged on a sphere, where with each edge a 
Hilbert space "Hi oc is associated (generalizations to other 
2D lattices are possible). Finally we assign an arbitrary 
orientation to each edge of the lattice. Now consider the 
Hamiltonian i? qc j = — Bp ~ Yl A; where the first (sec- 
ond) sum runs over all plaquettes p (vertices v) of the 
lattice. The operators B p and A v are mutually commut- 
ing projectors with support along the boundary of the 
plaquette p and on the edges adjacent to v, respectively. 
To define B p , consider the edges (ei, e 2 , e 3 , e 4 ) of p when 
traversing this plaquette counter-clockwise starting from 
some arbitrary vertex. The signature of the edge a is 
+1 if ei is pointing along the counter-clockwise sense, and 

— 1 if it is pointing in the opposite sense. Now consider a 
basis state |gi, <?2, ff3, 34) where gi £ G is associated with 
the system on edge e,. Then B p acts as follows: 

B p \gi,g2,93,94,) = S{g s 4 4 g s 3 3 g2 2 g s 1 1 )\gi,g2,g3,g4), (A2) 

where 5(e) = 1 and 8(k) = for every e 7^ k £ G (where 
e is the neutral element of G). To define A v consider 
the four edges (/1, /2, /3, A) incident on v. The signa- 
ture ti of fi is now +1 if fi is pointing towards v and 

— 1 if it is pointing away from v. Consider a basis state 
101) 9ii ff3, 94) where gi £ G is associated with the system 
on edge fi. If k £ G define 

U=-i (A3) 

Now define 

A,|si,ff2,ff3,ff4) = 7777 Y] \ k *9i, k *92,k*g3,k*g 4 )- 
|G| fcGG 

(A4) 

The ground state subspace 7W q d of i? q d consists of those 
\4>) which are joint +1 eigenvectors of all projectors B p 
and A v . To show that M q d is an M-space, define 

U p = 2B p -I 
V v (k)\gi,g 2 ,93,9i) = \k * gi,k * g 2 ,k * g 3 ,k * g 4 ), 

(A5) 



where k £ G. Then U p is a diagonal unitary operator 
and V v (k) is a permutation operator. It is immediate 
that the +1 eigenspaces of U p and B p coincide. Further- 
more the +1 eigenspace of A v coincides with the +1 joint 
eigenspace of the set A v := {V v (k) : k £ G}. To see this, 
note that (a) A v is a group and (b) A v is the averaging 
operator over this group. It follows that (cf. appendix 
[Bj) A v is the projector onto the +1 common eigenspace 
of A v . In other words the +1 eigenspace of A v coincides 
with the +1 joint eigenspace of A v . 

In conclusion, Ai q d coincides with the +1 common 
eigenspace of the operators U p and V v (k) where p ranges 
over all plaquettes, v ranges over all vertices and k over 
the entire group G. This shows that M q d is an M-space. 



3. Laughlin wavefunction at v — 1 

The Laughlin wave function for n particles and filling 
fraction v = 1 is 

n-l 

|^)= ]T e a ^\ ai )...\a n ) (A6) 

eti,...,ein = 

Here \ai) is a wavefunction for particle i given by 

(z\a % ) - -l=^ e -N 2 /2 (A7) 
V 7ra! 

where z is a complex number encoding the position of 
the particle in a two-dimensional plane. Furthermore e 
is the Levi-Civita completely antisymmetric tensor in n 
dimensions. Let Sij be the operator which swaps systems 
i and j. It is straightforward that each operator —Sij is 
unitary and monomial and that |^) is the unique +1 
common eigenvector of these operators. Thus is an 
M-state. 



4. Coherent probabilistic computations 

Consider a poly-size classical circuit composed of re- 
versible gates (e.g. Toffoli gates and NOT gates). The 
circuit acts on an n-bit input string where, say, the first 
k bits are uniformly random and the last n—k input bits 
are initialized in 0. Let {n x : x £ l^} denote the out- 
put probability distribution over the set of n-bit strings 
and define the state |tt) = J] ^JW^\x). Let C denote the 
natural translation of C into an n-qubit quantum cir- 
cuit. Furthermore let Xi and Zj denote the Pauli X and 
Z operators acting on qubit i and j, resp., and define 
Pi = CXi& and Dj = CZjC* . These operators are uni- 
tary and monomial (relative to the computational basis) 
and have \tt) as unique +1 common eigenstate. Thus \w) 
is an M-state. 
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5. Coset states of finite Abelian groups 

Consider a Hilbert space with basis {\g) : g € G} 
where basis vectors are labelled by the elements of a finite 
abelian group G. Let \g be the character of G canoni- 
cally associated to g and define the operators 

X(g): \ x )^\x + g); Z(g): \x) -> Xa(x)\x), (A8) 

for every g, x G G. Remark that these operators are 
unitary and monomial. Consider a subgroup H of G 
with generating set {h 1 , . . . , h n }. The dual group H is 
the set of all k satisfying Xh(k) = 1 for every h G H; 
consider a generating set {fc 1 , . . . , k" 1 }. Now define the 
"coset state" \H+x) = \x+h) where the sum is over all 
h G H. This state is the unique +1 common eigenvector 
of the operators X(h l ) and Z{k 3 ) where i ranges from 1 
to n and where j ranges from 1 to m. Thus such coset 
states are M-states. Remark that these states play an 
important role in the context of quantum algorithms viz. 
the Abelian hidden subgroup problem. 

6. W-states and Dicke states 

Consider a system of n qubits with computational basis 
B. Fix fee {1, . . . , n} and define \D k ) to be the (properly 
normalized) equal superposition over all basis states \x) 
for which the bit string x has hamming weight k. Note 
that ID 1 ) = \W) is the W state. Denote T k = a k ■ A®", 
where A = diag(l,a) and a = e~ . Further, let Si 
denote the SWAP gate acting on qubits i and i + The 
operators Tk and Si are unitary and monomial and have 
\D k ) as unique +1 common eigenstate. Thus \D k ) is an 
M-state. 



Appendix B: The projector p 

Let 7r = \Q\~ X Y2,U be the averaging operator over Q. 
We prove that tt — p. Since tt is defined as a sum over 
all elements in Q, one has Uir = tt for every U € Q. 
Consequently, 

^ = w\ £ Un = *• ( B1 ) 

11 ueg 

Moreover, tt = since for every U G Q one has G Q 
as Q is a group. Thus tt 2 = tt = tt> i.e. tt is an orthog- 
onal projector. Now let W denote the +1 eigenspace of 
tt. First note that the definition of tt immediately im- 
plies that Tr\ip) = for every G A4, showing that 
A4 C W. Furthermore, consider an arbitrary \ip) G W 
i.e. 7r|i/3) = |</?). Since Utt = tt for every U E Q. 'it follows 
that 

U\<p) = Uir\y) - tt\ V ) = \<p), (B2) 

showing that \<p) G A4. This proves that W C A'i and 
we conclude that W = A^. Thus p and 7r are orthogonal 



projectors on the same space, so that these operators are 
equal. 

To prove consider an arbitrary fixed U £ Q. Since 
Q is a group, the sum XVeG ^ s P rec i se ly the sum 
over all elements of Q. Using ([6]) it follows that Up = p. 

Appendix C: Proof of theorem [3] 

Problem 1 is shown to be NP-hard by relating it to 
3-SATISFIABILITY (3SAT). In the latter problem, the in- 
put is a Boolean function f(x) (where X - — - 3/1 * * * IS 
a bit string) written in conjunctive normal form with 3 
variables per clause (3-CNF). That is f(x) = fi(x) A 
■ ■ • A fk(x) where each clause /, is a disjunction of three 
literals; a literal is either a variable or its negation. An 
example of a clause is X3 V V X7. The 3SAT problem 
is to decide whether / is satisfiable i.e. if there exists 
an x such that f(x) = 1. To reduce 3SAT to Problem 
1, consider an n-qubit system and for each i define the 
unitary operator [7j which maps \x) to itself if fi(x) = 1 
and to —\x) otherwise. It is then easy to verify that the 
operators U\, ■ ■ ■ , Uk have a +1 common eigenvector if 
and only if / is satisfiable. Moreover every Ui is unitary 
and diagonal an acts on at most 3 qubits. 

As for Problems 2 and 3 we consider the following vari- 
ant of 3SAT. The input is again a Boolean function / in 
3-CNF form; in addition it is promised that there exists a 
unique x* such that f(x*) = 1. The problem is to deter- 
mine x* . This problem is known to be NP-hard (under 
randomized reductions) 25]. To reduce it to problem 2, 
consider the translations of the clauses fi into diagonal 
unitary operator Ui as above. Owing to the promise on 
/, it follows that the state :— \x*) is an M-state with 
stabilizer group (U\, ■ ■ ■ , Uk)- Sampling the distribution 
II is equivalent to determining x* . 

As for Problem 3, note that estimating the n expec- 
tation values (tl>\Zi \ip) in polynomial time with the pre- 
scribed accuracy allows one to determine x* in polyno- 
mial time as well. 



Appendix D: Proof of lemma [3] 

Clearly, the procedure in the lemma can only output 
elements lying within O x . Note that every \y) in this 
orbit occurs with probability 

Prob(|y» = \{P G V : \y) = P\x)}/\Q\. (Dl) 

Let Q consist of all P G V satisfying P\x) = \x). Then 
Q is easily shown to be a subgroup of Q. Further, let 
Py G Q be a permutation satisfying P y \x) = \y) and 
consider the set P y Q consisting of all products P y P with 
P G Q (formally, this is the left coset of Q defined by 
P y ). We now claim that 

P V Q = {P G V : P\x) = \y)}. (D2) 
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The inclusion C holds trivially. To prove the converse 
inclusion, note that every P satisfying P\x) — \y) can be 
written as P = PylP^P], where Py X P £ Q. This shows 
(|D2l) . Therefore Probfly)) is equal to \P y Q\/\V\. Note 
however that by definition the set P y Q has the same car- 
dinality as Q for every y. Thus Prob(|y)) is independent 
of y so that this probability distribution is uniform, as 
desired. 



Appendix E: Pauli operators 

For every a, a', x £ Z 2 , the following well-known iden- 
tities are easily verified: 



X(a)\x) 


= \x + 


a) 


(El) 


Z(a)\x) 


= (-i; 


xTa \x) 


(E2) 


X(a)X(a') 


= X(a 


+ o') 


(E3) 


Z(a)Z(a') 


= Z(a 


+ a') 


(E4) 


X{a)Z{a!) 


= (-i; 


aTa 'Z(a')X(a) 


(E5) 



If a = i k X(s)Z(t) (where k £ {0,1,2,3} and s,t £ 1%) 
is a Pauli operator we will call (k, s, t) the label of a. 

Proof of lemma^ we will need take the following tech- 
nicality in account. Remark that a 2 — al for every 
Pauli operator a, where 7 = ±1 . Now suppose that 



there exists a £ Q where 7 is not equal to 1. Then, since 
7/ = a 2 £ Q one would have 7/^) = \ip) so that \ip) = 0, 
leading a contradiction. Thus we conclude that a 2 = I 
for every a £ Q . 

Since the o~i mutually commute and square to the 
identity, every operator in Q can be parameterized as 
a(a) := cr" 1 . . . <r"", where a = (ai, . . . , a n ) £ ZJ. Using 
(|EliE5|) it follows that 

a(a)(xX(J2a t s t )Z{J2a l t l ) (E6) 

This implies that a (a) £ T> if and only if a solves the 
linear equation ^a^s* = 0- Let {a 1 , . . . ,a k } denote a 
basis of solutions; such a basis can be computed effi- 
ciently. Thus <j (a) £ T> if and only if a — ^2yja J for 
some yj £ Z2. Furthermore, using the definition of o~(a) 
it is straightforward to show that 

o-(El/j ffli ) = ^(a 1 )^ .. .a{a k ) v K (E7) 

This shows that the operators Dj :— cr(a- ? ') form a gen- 
erating set of T>. Since o~(a J ) is a product of at most n 
Pauli products, the label of each a {a?) can be computed 
in poly(n) time. Finally, since every element of Q squares 
to the identity it follows that Dj = (— 1)"^' Z(d?) for some 
uj £ Z 2 and some d? £ Z£. □ 
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